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Abstract 

The energy-weighted sum rule for an electric dipole transition operator of a Schiff type differs 
from the Thomas-Reiche-Kuhn sum rule by several corrective terms which depend on the number 
of system components, J\f. For illustration the formalism was applied to the case of Na clusters. 
One concludes that the RPA results for Na clusters obey the modified TRK sum rule. 
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I. INTRODUCTION 

The Thomas- Reiche-Kuhn (TRK) sum rule [l^l has been widely used in various con- 
texts of electron excitations in atoms, molecules, and solids. Indeed, the relation to the 
photoabsorbtion oscillator strength makes it quite useful in interpreting various collec- 
tive features of the dipole states {]- 1(3]. 

The many body properties like collective excitations in atoms, molecules or clusters are 
studied using different collective models. Thus, the random phase approximation (RPA) 
formalism has been extensively used and excellent results for the photoabsorbtion cross 
section spectra, especially in the atomic clusters domain [8 

In a previous publication [lo|], some of the many body features of the small and medium 
sodium clusters were studied within the RPA approach using the projected spherical single- 
particle basis defined in Ref . jll| . The RPA wave functions were used to treat the dipole 
transitions which led to the photoabsorbtion cross section spectra. Also, the system static 
electric polarizability was analytically determined. 

The salient feature of the RPA approach with spherical single particle basis consists of 
that it satisfies the TRK sum rule. This is however true for the electric dipole moment, 
which is not the case in the above quoted paper where, indeed, a modified dipole operator, 



10| were obtained. 



similar with the Schiff dipole moment 
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16] is used. Using the standard dipole operator 



our formalism is not able to describe the dipole excitations beyond an energy of about 3eV, 
because one cannot account for the transitions between states characterized by AiV = 3. 

The question which, however, should be answered is whether the use of the Schiff momen- 
tum is imposed by a physical constraint or by a technical difficulty We attempt to provide 
the answer below. The dipole moment manifests itself, for example, in the interaction of 
atoms with external electric field which results in inducing a Stark energy shift. A two body 
dipole interaction of dipole-dipole type violates the conservation of the center of mass mo- 
mentum. Treating the interacting system of electrons moving in the electrostatic potential 
created by the ionic core, and an external electric field in the reference frame where the total 
momentum of the electron system is conserved one obtains two important results: 1) The 
total electric field acting on the ionic core is equal to zero, this result being known under 
the name of the Schiff theorem. 2) The interaction of electrons with the core is described 
by the electrostatic potential, generated by the ionic core, and a screening term which is 
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proportional to the gradient of the first. The total potential energy can be written as a 
series of the average values < r n >. The first non- vanishing P, T-odd term of this expansion 
is a mixture of the octupole and the Schiff momenta. Therefore, the correction to the dipole 
appear to be a recoil term which assures the restoration of the conservation law for the center 
of mass linear momentum. In conclusion, the correction to the dipole operator accounts for 



the screening effect, caused by the electronic shells, on the motion of the valence e 



ectrons 



10| uses a 



in the mean field determined by the ionic charge. Moreover the approach of Ref. 
projected spherical single particle basis which allows for an unified description of spherical 
and deformed clusters. 

Here we address the question whether the specific features of our approach require a 
modification of the TRK sum rule. The result of our investigation is that a new sum rule 
holds when one uses the Schiff-like dipole momentum in order to describe the photoabsorbtion 
spectra in Na clusters. 

In what follows we shall show how did we arrive at this interesting result. First we provide 
few basic information. 



II. THE ENERGY WEIGHTED SUM RULE 

Within the RPA formalism, for any Hermitian operator M, the following sum rule holds 
^2(E n - Eq) |(0|M|l n )| 2 = ~<0|[[M,#],M]|0), (2.1) 

n 

where E n are the RPA energies associated to the many body Hamiltonian H. Here the state 
|0) is the RPA phonon vacuum, while |l n ) denotes the single phonon state |l n ) = C^|0). 

This result is known under the name of Thouless's theorem [nj]. The summation in 
Eq. (12 .ip suggests that the sum rule may give information about the global properties 
of the many body system. Indeed, by using the the lowest k-th energy weighted moment 
of the strength function one obtains simple expressions for the mean energy (or energy 
centroid) and the variance of the strength function for the linear response of the system 
to an excitation operator. Lower and upper bounds for the mentioned quantities can be 
analytically obtained. Also the static polarizability of a cluster is proportional to the sum 
rule S-2- 

We are interested in those properties of atomic clusters which are determined by the 
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motion of the valence electrons. These move in a mean field, determined by the ionic core, 
and interact among themselves through a Coulomb force. The residual two-body interaction 
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191 ] can be expanded in multipole series from which only the dipole term is relevant and 



therefore considered. The phonon operator is defined as: 

cftM = £ [x; h (4c h ) lfl - 1^(4^] , (2.2) 

ph 

where X£ h and Y™ h are the n-th order solution of the RPA equations, cj, and denote 
the creation operators for particles and holes, respectively. The reduced probability for the 
dipole transition |0) — > can be written in terms of the RPA phonon amplitudes and the 
ph matrix elements of the transition operator M : 

B(El,0+^l-) = \(0\\M(El)\\l-)\ 2 . (2.3) 



Instead of the usual transition dipole operator, a Schiff-like moment operator 
was used. 



3, HQ 



M{El) = e ( 1 - -— ) f. (2.4) 



3r_ 2 ' 
5r 2 

Here r s is the Wigner-Seitz radius which, for Na clusters has the value of 3.93 a.u.. The 
corrective component, involved in the dipole operator, may relate particle and hole states 
characterized by AiV = 3, which results in modifying the strength distribution among 
the RPA states. Such an effect would be however obtained even for the dipole transition 
operator, if the mean field potential for the single particle motion involves higher powers of 
the radial coordinate. 

Reckoning the double commutator from Eq. f)2.ip . corresponding to the transition oper- 
ator (12. 4p one obtains: 



qt.2 2 

J2(E n ~ E )\(0\\M(El)\\l n )\ 2 = — 

n e 



X 



^-^(oiiE^i ) + ^i( iiE^i°) 



(2.5) 



The terms correcting the standard TRK sum rule are the expected values of the radius 
powers r 2 and r 4 , in the RPA ground state. Obviously, these corrective terms induce 



[1] Throughout this paper the Rose's convention for the reduced matrix elements are used. 
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an additional M dependence for the energy-weighted sum of the reduced dipole transition 
probabilities. In what follows we shall use the notation EWS for the left side and S(N) for 
the right hand side of Eq. ( 12. 5p . EWS can be directly calculated using the RPA output data, 
like energies and transition probabilities [l0|. If S(Af) has a model independent expression 
involving only universal constants, and the equality EWS = S(J\f) holds, one says that a sum 
rule is valid for the system under consideration. Since EWS is calculated by complicated 
many body formalisms, obeying the sum rule equation is a serious test for the proposed 
description. The terms of S(Af) were alternatively evaluated through two distinct methods: 



A. The boson expansion method 



The terms involved in S (Af) can be expressed in terms of particle-particle (pp) and hole- 
hole (hh) transition matrix elements. Indeed, the ph transition components give vanishing 
contributions when they are averaged with the RPA ground state. Therefore, the needed 
one-body operator f m can be written as: 



Af 



E r * = E^^K^+E^^i^- (2.6) 

o=l p h 

The fermion density operators c p c p and c\ch can be expressed in terms of the RPA phonon 
operators C\{C n ): 

n n 

where the coefficients a™ and b 1 ^ have the expressions: 

^ = (0|[[Ct, c t Cp ],C n ]|0), bl = (0\[[Clclc h },C n }\0). 

In this way the modified dipole sum rule becomes: 

, 2 



J2(E n - E ) \(0\\M(El)\\l n }\ 2m 



(2.7) 



(2.8) 



(2.9) 



■ph ) 



n,p,h 
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~2 r (P h ) + W 



where the following notation has been used: 

r k (ph) = (p\r k \p) - (h\r k \h), k = 2,4. (2.10) 

As we have already mentioned, a specific ingredient of the approach from Ref. 10[ is the 
use of the projected spherical single particle basis in the RPA formalism. 



B. Electron density approach 



Static electric polarizabilities results based on the calculations for the number of spilled 



out electrons 



10], agree quite well with the corresponding experimental data. These spilled 



out electrons produce a screening effect against external fields which results in changing the 
classical result for the polarizability. The basic assumption in accounting the spilled out 
electrons is the fact that the electron density is not going sharply to zero at the cluster 
surface, but is gradually decreasing and moreover extends significantly beyond the jellium 
edge. The same argument can be brought for the correction terms of the S(Af) containing 
averages of the radius powers with the RPA vacuum state. The electronic density has a 
constant central part, enfolded by a diffuse region, of a width equal to a, where the electron 
density tends smoothly to zero. Guided on some formal parallelism between the behaviors 
of atomic clusters and nuclear systems 20, |2l|, the average of r m with the RPA vacuum 



state is approximated by folding r m with a localization probability density with spherical 
symmetry, of Fermi distribution type. 

/ r — n \ 

Po 



p(r) 



1 + exp 



(2.11) 



Here R = r s M 1 ^ is the radius of the cluster with Af atoms, a is a parameter defining 
the thickness of the diffusion region. For nuclear systems, a is a constant related to the 
diffuseness parameter d by a simple relation: d ~ 4.4a. 



For atomic c 
and Sorbello 



usters such a density function was used, in a different context, by Snider 



3 if ] - A modified density function, having the denominator at a power 7 / 1, 



which yields a density profile with an asymmetry around the inflexion point, was used in 



Ref. 



19|. 



Having the density, momenta < r m > may be written in the form of a power series in the 



variable J| 



22|: 



(r m ) 



R r 



m + 3 
m = 2,4. 



7r 2 /a\ 2 . . 

T \r) m ( m + 5 ) + --- 



(2.12) 
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In this way the term S(Af) can be written in the following way 



9h 2 e 2 
2m 
7i 2 a 2 /594 



M-^N 2 l* + ^N^+ 



5r? V 35 



s 



A/" 2/3 - 14 



2^2 



9h z e 
2m 



-J 7 (AO- (2.13) 



One notices that besides the number of atoms dependency, this expression involves only 
universal constants, which makes Eq. (12.13!) to be, indeed, a real sum rule. 

Contrary to the case of nuclear systems, where the thickness of the diffusion region 
is approximately the same for all nuclei, for atomic clusters the parameter a is expected 
to have a M dependency due to the long range character of the two body interaction. 
This dependence was determined by interpolating the values of a satisfying the equation 
EWS = S(J\f), for 8 < M < 40. The solutions for a were interpolated by the function of 
J\f: 

a(J\f) = -0.975 - 0.0lW 1/3 + 0.361A^ 2/3 . (2.14) 

Inserting the expression (I2.14p of a(AT) in Eq. (12. 13[) . one obtains the final expression for the 
sum S(Af). 

The sum rule is always a serious consistency test for any model calculation and in partic- 



ular for the RPA formalism of Ref. jlO| which uses a projected spherical single particle basis, 
appropriate for an unified description of spherical and deformed clusters, and a Schiff-like 
momentum as a dipole transition operator. 

Since the double commutator appearing in the sum rule (2.1) is averaged with the ground 
state, it results that the sum rule reflects the structure of the ground state. 

We would like to mention that Eq.f l2.ll) is valid also when the sates involved are exact 



321 ] . Note that despite the fact 



eigenstates while energies are the exact eigenvalues of H 
RPA is an approximative approach the TRK sum rule is exactly obeyed. Due to this feature 
it is a legitimate question which are the approximations which preserve the sum rule and 
which are those which violate it. Of course this question cannot be completely answered. 
Despite of the fact that this issue is more suited for a revue paper rather than for a short 
letter, here we give some brief comments on this matter. 

If we replace the RPA dipole states by the particle-hole dipole states the TRK sum 
rule is still satisfied. In nuclear physics the N-Z sum rule for the Gamow- Teller transition 
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operator is also satisfied in the framework of BCS approximation. In this case the proton- 
neutron quasiparticle RPA dipole states are replaced by two proton-neutron quasiparticle 
dipole states. Aiming at satisfying the Pauli principle, the pnQRPA has been improved by a 



renormalization procedure 
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34| , while the transition operator was left unmodified. With 



these new ingredients the N — Z sum rule is violated by about 25-30%. Another approach for 
this transition process was formulated by one of the present authors (A.A.R.) [35j, by taking 
for the transition operator a first order boson expansion expression while the states involved 
for the mother and daughter nuclei were still described within the pnQRPA formalism. 
Again the deviation from the sum rule is about 20-25%. It is interesting to remark that 
the first procedure underestimates the sum rule while the second one overestimate it. This 
suggested us to use for the GT transition operator a boson expansion expression in terms 
of the pnQRPA renormalized boson. The states assigned for the nuclei participating to 
transitions are described by the renormalized pnQRPA formalism. In this way the deviation 
of the sum rule from the N — Z value was diminished up to less than 10%. 

An inconsistency which might induce a deviation from the TRK sum will be discussee 
a bit latter, when the numerical results for the boson expansion method applied to the 
expected values of r 2 and r 4 , are presented. 

The structure of the mean field may be a source of inconsistency. Indeed, if the mean field 
is of a harmonic oscillator type the dipole operator is not able to connect states characterized 
by AN = 3 and consequently the volume like collective excitation would be missing. 

An important feature related to the sum rule is the completeness of the states defined 
within the chosen approximation. For example, the RPA calculations are performed within 
a restricted single particle and hole space which determine the unity resolution for the RPA 
space( 1 = J^fc |lfc)(lfcl) which, in its turn, yields a certain value of the sum rule. If we 
enlarge the single particle space, the new RPA space determines a new unity resolution 
and consequently a new value for the sum rule. In general, the sum rule depends on the 
dimension of the single particle space used by the RPA approach. In principle the sum rule 
is an increasing function of the dimension of the particle-hole configurations. How fast the 
sum rule converges when this dimension is increased depends of course on the single particle 
basis provided by the employed mean field as well as on the transition operator. 

The exchange term involved in the electron two body interaction yields a screening effect 
for the system response to the action of an exciting external field. However such a screening 
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effect is not found in the transition dipole operator. As a matter of fact this is the incon- 
sistency we wanted to remove when we used the Schiff-like dipole operator. Moreover, such 
an operator links the AiV = 3 states and, therefore, yields a non-vanishing strength for the 
volume dipole excitation. 

In our formalism the RPA procedure uses a projected spherical single particle basis and 
a Schiff-like dipole operator. Due to these new ingredients we expect a new sum rule. As we 
have seen before, the corresponding EWS has a simple expression as a function of M and 
consequently satisfies the commonly used definition for a sum rule. The consistency of our 
calculations is tested by checking numerically whether the equation 12.101 is satisfied. 

III. NUMERICAL APPLICATION 

The EWS was calculated with the RPA energies and matrix elements for the transition 
operator. S(Af) was alternatively calculated with the expressions (9) and (13), respectively. 
The two sets of results corresponding to the mentioned options for S(Af), are plotted in 
Figs. 1 and 2, respectively. 

The RPA calculation of S(Af) does not yield an explicit Af dependence although such 
a dependence is involved in the many body formalism by means of the Fermi energy, the 
oscillator length, the matrix elements in the projected single particle basis and the space 
of the single particle states used in the RPA calculation. Note that Fig. 1 shows a good 
agreement for medium clusters. The deviation is significant for large as well as for small 
clusters although both curves exhibit similar pattern concerning the oscillating behavior. 
The discrepancies may indicate that higher order boson expansion terms are necessary in 
order to improve the agreement. However, these terms would bring a certain inconsistency 
to the formalism since EWS is evaluated within the RPA approach. Moreover, in order that 
Eq.(9) plays the role of a sum rule it is necessary that S(J\f) exhibits a model independent 
expression. 

In this context it is worth mentioning that the second procedure described above, yields 
indeed, a model independent expression for S(Af). The results for this situation are shown 
in Fig. 2 where a very good agreement between S(J\f) and EWS is shown. 

Actually this is the main result of our investigation. The sum rule for the modified dipole 
operator involves S(Af) given by Eq.(2.13) with the parameter a from (2.14)- By contrast 
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FIG. 1: The calculated EWS (open circles) and S{M) given by the RPA approach (black triangles), 
versus J\f, the number of cluster's components. 
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FIG. 2: The same as in Fig. 1 but S{J\f) given by Eq. (13) and represented by a solid line. 
Theoretical result a given by black circles. 

to the TRK sum rule, the new one is not linear in N '. 

Of course one may ask how the sum rule is mirrored by the experimental results. Direct 
data for the momenta < r 2 > and < r 4 > are not yet available. Under these circumstances 
we have to see whether the experimental data for the EWS agree with the calculated 
values. Actually, the experimental value for EWS can be extracted from the experimental 
photoabsorbtion cross sections. Indeed, interpolating the discrete values of the experimental 
photoabsorbtion cross sections, given as function of the excitation energy, by a smooth curve 
and integrating the result with respect to the energy in the interval [0, oo), one obtains the 
area A(J\f). The experimental EWS is proportional to A(J\f). We assume that the M 
depending proportionality factor is the same as for the theoretical EWS, i.e. J-(Af). This 
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assumption is grounded on the fact that F{M) is a model independent quantity and moreover 
assures similar normalization for the total cross section as in the schematic calculation. Thus, 
the experimental EWS is defined by: 



[EWS] 



Exp. 



(3.1) 



The quantity A(Af), extracted from the data of Refs. 28_|, |29|], varies between 0.582 (eV)A 2 
for M = 9 and 0.387 (eV)A 2 , for N = 19. Here the constant factor is Q = 70.159[e 2 ]. This 
value was obtained by equating the right hand side of the above equation to the calculated 
value for EWS for the case of Nau, were the best agreement between the calculated and 
experimental polarizabilities was obtained. Thus, the constant Q yields a normalization of 
EWS which accounts for the "missing sum rule", noticed experimentally 30]. 

Actually, this normalization of the integrated cross section might be used for certifying 
the collective character of a dipole excitation, according to its contribution to the sum rule. 
Thus the collective surface dipole excitations exhaust between 70% to 100% from the sum 
rule. In the case of small and medium clusters the remainder sum rule is attributed either 



to a volume-like excitation [26[ or to single electron excitations 29]. 




FIG. 3: The quantities [-EW5] Exp. (crosses) and [EWS]th. (open circles) defined by Eqs. (15) and 
(16) respectively, are given as function of the number of cluster's components, M. 



In Fig. 3, the quantity [EWS] Exp is compared with the theoretical weighted sum EWS, 
calculated within the RPA approach, which uses a projected spherical single particle basis 



[EWS] 



Th. 



J2(E n - e )b(o+ -> 1-; 



(3.2) 



RPA 
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The agreement between the two EWS's, shown in Fig. 3, is a guarantee that [EWS] E is, 
indeed, close to the S(M). Note that in Fig.3 we didn't consider clusters with 21 < Af < 40 
although for some of them experimental data are available. The reason is that for these 
data only few cross sections are given and therefore it is not possible to calculate the area 
A introduced above. 



IV. CONCLUSIONS 

Finally, we want to mention that the octupole correction to the dipole transition oper- 
ator was also used in connection with the description of the electric dipole transitions in 



nuclear systems. Thus, in Ref . [23J it is pointed out that adding the octupole correction, the 
agreement with experimental data concerning the El transitions is substantially improved. 
In Ref. 24| a similar effect is obtained by modifying the many body wave function due to the 
octupole interaction which is considered in addition to an isovector-dipole interaction. The 
new components are connected by the standard dipole operator and consequently modifies 
the El transition rates. Of course, an energy-weighted sum rule associated to the dipole 
transition operator holds also for nuclear systems. Obviously, changing the transition op- 



erator, as it happened in Ref.j23|, the corresponding sum rule which should be valid is the 
one obtained in the present paper. 

Another interesting example is the N — Z sum rule, which holds for the single beta 
transition. This rule says that for a single beta decaying nucleus, the difference between 
the f}~ and j3 + strengths should be equal to 3(iV — Z). One remarks that this is in fact 
the nuclear physics counterpart of the TRK sum rule, which may be alternatively formu- 
lated as expressing the equality of the difference between the photoabsorbtion and emission 
dipole transition probabilities and the number of delocalized electrons. The particle-hole 
RPA approach satisfies exactly the TRK sum rule. The N-Z sum rule is true, for example, 
for Gamow- Teller (GT) dipole transitions and is exactly satisfied within the proton-neutron 
quasiparticle RPA approach. Extensions of the microscopic formalisms to the double beta 
decay 2i/(3f3, showed that in order to describe the transition rates, it is necessary to im- 
prove the wave functions of the mother nucleus as well as the GT dipole states, by adding 
anharmonic effects. However, these corrections violate drastically the N — Z sum rule. In 
the spirit of the present paper, we open the question whether the Gamow- Teller proton- 
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neutron interaction could be extended by adding an octupole component such that to the 
new proton-neutron interaction a modified N — Z sum rule corresponds. This would make 
the inclusion of anharmonic effects which, as a matter of fact, violates the Pauli principle, 
unnecessary. It is an open question whether a Schiff like correction of the Gamow- Teller 
transition operator is necessary due to some specific conservation law, or just due to the 
necessity of improving the existent descriptions. 

The final conclusion is that the Schiff-like dipole moment used for the RPA description 
of the photoabsorbtion cross section spectrum, satisfies an extended TRK sum rule. The 
saturation of the extended sum rule is a positive test for the single particle basis as well as 
for the dimension of the dipole ph space involved in the RPA description. The TRK sum 
rule is of a general interest, being applicable also for other many body systems correlated by 
a Schiff-like two body interaction. The usefulness of sum rules in exploring the many body 
properties mirrored by the multipole electric, or magnetic transitions have been stressed by 
many authors 2514271] . Here we showed that a sum rule may hold also for a multipole mixed 
transition operator. 
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